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We study interacting bosons in optical lattices in the weak-tunneling regime in systems that exhibit 
the coexistence of Mott-insulating and condensed phases. We discuss the nature of the condensed 
ground state in this regime and the validity of the mean-field treatment thereof. We suggest two ex- 
perimental signatures of condensate order in the system. (1) We analyze the hyperfine configuration 
of the system and propose a set of experimental parameters for observing radio-frequency spectra 
that would demonstrate the existence of the condensed phase between Mott-insulating phases. We 
derive the structure of the signal from the condensate in a typical trapped system, taking into ac- 
count Goldstone excitations, and discuss its evolution as a function of temperature. (2) We study 
matter-wave interference patterns displayed by the system upon release from all confining potentials. 
We show that as the density profiles evolve very differently for the Mott-insulating phase and the 
condensed phase, they can be distinguished from one another when the two phases coexist. 

PACS numbers: 37.10. Jk,32.30.Bv,37.25.-|-k 



I. INTRODUCTION 

Ultracold bosons in optical lattices allow extensive 
studies of quantum phases of condensed matter in con- 
trolled environments [l|, 0] . When the tunneling energy, 
relative to the interaction energy, decreases across a crit- 
ical value, the system undergoes a transition from the 
superfluid phase to the Mott-insulating phase 0]. This 
transition has been observed in experiments by imaging 
matter wave interference patterns of the system In 
realistic systems, a confining trap renders the density of 
bosons non-uniform and, in a sufficiently deep optical 
lattice (where the lattice depth is over 20 times the re- 
coil energy), this inhomogeneous system is predicted to 
have a multiple-layer structure in which Mott-insulator 
layers having different occupation number are separated 
by relatively thin condensate layers Q. The study 
of such inter-layer condensates is important for a vari- 
ety of reasons. From the perspective of quantum in- 
formation schemes 0, i, H, robust, high-fidelity Mott- 
insulator states are required and the presence of a con- 
densate is in fact a hindrance. The condensate is also of 
interest in and of itself; its ground-state properties and 
excitation spectrum have been studied in various ways 
0, M, E M, Q M, M, [13 and show vast differences 
from more commonly encountered condensates in free 
space or in the strong-tunneling limit in shallow opti- 
cal lattices. Finally, the interlayer structure is of interest 
as a model realization of spatially coexisting quantum 
phases separated by critical crossover regimes. 

The purpose of this work is twofold: to provide a com- 
prehensive characterization of condensates in the vicinity 
of Mott-insulators and to propose settings and measure- 
ments which would pinpoint direct experimental signa- 
tures of the presence of condensate interlayers in the in- 
homogeneous Mott-insulator-condensate shell structure. 
We begin by deriving key properties of the interlayer 
condensate, such as off-diagonal elements of the single- 



particle density matrix, number fluctuations, and con- 
densate fraction, in the deep-lattice (weak-tunneling) 
regime. We describe the condensate ground state with 
a fixed-number wave function, which displays entangle- 
ment and off-diagonal long-range order, and conserves 
the total number of particles as any exact eigenstate 
of the Bose-Hubbard Hamiltonian should do. We com- 
pare this fixed-number entangled wave function with a 
commonly-used mean-field wave function. By calculating 
the condensate properties in each description, we iden- 
tify the conditions under which the mean-field approxi- 
mation is valid. These studies shed light on the range 
of validity of mean-field descriptions of the condensate 
in the deep-lattice regime. We analyze the nature of the 
excitations of these condensates, including the spectrum 
of Goldstone mode excitations, that would be relevant in 
the oft-encountered experimental situation of two species 
of bosons. 

Recent experiments using hyperfine transitions driven 
by external radio- frequ ency (rf) drives have conclusively 
shown the existence [3, US and formation of a Mott- 
insulating shell structure in systems where the conden- 
sate interlayers are predicted to be much thinner than the 
Mott-insulating shells. Currently, however, no strong ex- 
perimental evidence unequivocally verifies the existence 
or properties of the putative condensate layers between 
the Mott-insulator layers [l^. In this work, we propose 
a specific set of parameters for an rf experiment in which 
the system is prepared and resolved as in Ref. [l^. We 
find that the presence of the condensate interlayer can 
lead to a two peak structure, in contrast to a single peak 
structure associated with each adjoining Mott-insulator 
phase. We discuss how this signature rf spectrum pro- 
file ought to be robust against Goldstonc-mode pertur- 
bations associated with spontaneous symmetry breaking 
in the condensate and against low-temperature effects, 
and would be sensitive to the destruction of condensate 
order. This sensitivity to condensate order should allow 
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experiments to distinguish whether the incommensurate- 
density layers between Mott-insulator shells are in the 
superfluid or normal fluid phase. 

Another commonly used probe which has provided a 
wealth of valuable information on the Mott-insulator- 
superfluid transition is matter wave interference 
Here, we tailor the analysis of matter-wave interference 
to the specific case of the interlayer setting. We make 
transparent the contribution of various terms in the fixed- 
number entangled wave function to interference patterns 
that do not occur in the Mott-insulator state. While 
discerning the condensate amidst Mott-insulator phases 
through interference experiments can be a challenging 
task, our analysis provides methods to do so. 

The paper is outlined as follows. In Sec. |TT1 we review 
the Bose-Hubbard model and its parameters and detail 
two descriptions of the condensate ground state in the 
deep-lattice regime. We review how results for a uniform 
system can be applied to an inhomogeneous system using 
the local density approximation. In Sees. IIIII and llVi we 
propose feasible settings for experimental confirmation 
of the condensate layers in the inhomogeneous system at 
low temperatures using rf spectroscopy and matter-wave 
interference, respectively. In Sec. |Vl we summarize our 
results. 



II. CONDENSED BOSONS IN DEEP OPTICAL 
LATTICES 

A system of bosons in a deep optical lattice is well 
described by the Bose-Hubbard Hamiltonian 0, [2l[ : 



Hbh = ^ + ^1"' 



E 



^hi{n^ - 1) - (/i - Vi)hi 



(1) 



Here a,; and are the boson creation and annihilation 

operators on the ith lattice site and fii = a\ai is the 
number operator on the ith lattice site. ^ i denotes 
the sum over all single sites and ^ ^^j) denotes the sum 
over all nearest-neighbor pairs of sites. ^ is the chemi- 
cal potential, which is determined by the condition that 
the quantity '^i{nij is equal to the total number of 
bosons in the system. Vi is the value of the external 
confining potential on site i. For a homogeneous sys- 
tem, Vi is set to zero for convenience. J is the tunneling 
strength, which quantifies the ability of a boson to hop 
from one site to its adjacent sites. Assuming the number 
of particles per site is not too large {{fii) w 0(1)) 
J = - / d^Yw*{r - Ti) [-h^V'^/2m + Viat(r)] w{y - r^), 
where Viat(r) is the lattice potential, w{t — r.^) and 
w{v — Yi) are the single-particle Wannier wave func- 
tions of the lowest Bloeh band localized to the nearest- 
neighbor sites i and j, and m is a single boson's mass. 
U is the interaction between two bosons on a single site. 



U = {A-Kh^as/rn) J (Pr [^(r)!^, where is the scattering 
length of bosons, which is required to be small compared 
to the lattice spacing for the validity of this model. 

In experiments that do not use a Feshbach resonance 
to manipulate the interaction between bosons, the tun- 
able parameters are the lattice depth, Vq, and the lat- 
tice spacing, ds (half of the lattice laser's wavelength 
A), which in deep lattice are related to J and U by 
J = {4/^)ERiVo/Eiif^^eM-^VWE^] and U = 

?> / TTkiiasER{VQ / Ej^)'^/ , where = 27r/A and the re- 
coil energy Er = h^k\/{2m) 

In Sees. IIA-IIC, we discuss the condensate state in 
deep optical lattices. We discuss ground-state proper- 
ties of the condensate wave function which, as a result 
of its vicinity to Mott-insulator states, is distinctly dif- 
ferent from the condensate in the large-tunneling limit. 
We begin by constructing a ground state for a homo- 
geneous condensate in a fixed-number truncated basis, 
composed of superpositions of number states associated 
with the two neighboring Mott-insulating states. Using 
this ground state, we find expressions for the boson den- 
sity and the number fiuctuations on each site, calculate 
the single-particle density matrix, and hence obtain the 
condensate fraction. We find that the maximum conden- 
sate fraction in the deep-lattice regime is significantly 
smaller than one. We then consider a commonly em- 
ployed mean-field approximation to this state and show 
that the single-site number fluctuations, density, and con- 
densate fraction of this state are identical to those in the 
fixed-number state. We show that in describing the weak- 
tunneling condensate, using a truncated basis of occupa- 
tion number is valid, i.e., that the error induced by trun- 
cation is of order 0{J/U). We conclude with a discussion 
of the inhomogeneous system where Mott-insulator and 
condensate phases coexist. 



A. Fixed-number condensate state 

In the deep lattice regime(J/C/ ^ 1), the first term 
on the right-hand side of Eq. ^ can be treated per- 
turbatively to calculate physical quantities of interest, 
such as the single-site number fluctuation and the con- 
densate fraction, in orders of J/U. The unperturbed ba- 
sis is composed of all possible states of the form Yl ihh) 
with ^ rii equal to the total number of particles TV [23| ; 

is defined as (ni!)~^/^(a|)"' |vac). For convenience 
we use the symbol to denote a commensurate state 
that has each site occupied by n particles (nil"^))- The 
single-site number fiuctuation, denoted by Anf, is de- 
fined as Anf = ^n^) — (n^)^; the condensate fraction, 
denoted by /c, can be obtained from the single-particle 
density matrix of the system in the thermodynamic limit 
{N, M oo with N/M finite) [H]. Elements of the den- 
sity matrix are defined as (a|dj); the condensate fraction 
is the ratio of the largest eigenvalue of the density matrix 
to the total number of particles. 
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Given a system of M sites, if A'' is a multiple of M 
{N = nM, where n is an integer), the ground state of 
the system is A \ \n)) + 0{J/U), where the normalization 
constant A is equal to 1 up to order in J/U. For the 
state, Anf and fc are both zero. 

If = nM + All, where n and Mi arc integers and 
< All < M, the ground state of the system is 

I*) = E („+^;)!/./2 « ■ ■ • < II")) + 

_ (2) 

where {77} — {ryi, 772, ■ • ■ , VMi} ^ ^^t of distinct integers 
chosen from {1, 2, . . . , M}, and ^ {^j denotes the sum 
over all combinatorial configurations. The leading term 
shown in \'^) lies in a sub-space spanned by all possible 
product states of Mi single-site states with occupation 
number n + 1 and M — Mi singlc-sitc states with oc- 
cupation number n; in the small tunneling regime, this 
truncation proves to be a sufficient approximation. The 
coefficients C|^i, can be obtained by solving the Bose- 
Hubbard Hamiltonian, or cquivalcntly obtained by min- 
imizing the total energy under the normalization con- 

straint = 1- 

The number fluctuations on a single site, the density 
matrix elements, and the condensate fraction in the ther- 
modynamic limit can be readily calculated for this state 
(Appendix A). We find: 

An^ = (Afi/A/)(1 - Ml/A/), 

(a|a,;) = n + Mi/M, 

{a\aj^i) = (n+l)An^, 

fc = {a\dj^i)/{a\a,). (3) 

Thus, for the fixed-number state, the condensate fraction 
is the ratio of the off-diagonal elements of the density ma- 
trix to the diagonal one; the former is proportional to the 
number fluctuations and the latter is just the boson den- 
sity. For the fixed-number state, as stated in Eq.([3]), the 
off-diagonal elements are constant, reflecting a complete 
correlation within the entire system. 

If Ml = or Ml ~ M, An? and fc both vanish and 
the system is in the commensurate-filling Mott-insulator 
state \\n)) or 1)). If < Mi/M < 1, An^ and fc 

are both nonzero and the system is hence a condensate 
with a single-particle state macroscopically occupied; 
this state is the zero quasimomentum state (Appendix 
A). The condensate fraction has its maximum value 
{n + l)/(Vn + 1 + y/nf when Mi/M = yjn{n + 1) - n 
(the difference between n and the geometric mean of n 
and n + \). Compared to 100% condensation of the su- 
perfluid state in the shallow lattice limit (represented by 
\^i=ia\]^ \vac)), the condensate fraction in the deep 
lattice regime is less than 35% for all n> Q. 



B. Decoupled-site mean-field ground state 



In comparison with the fixed-number state above, 
we now investigate properties of a commonly employed 
mean-field ground state. The mean-field approach in- 
troduces a scalar field to represent the average effect of 
neighboring sites and leads to a site-decoupled Hamil- 
tonian but one which does not conserve particle num- 
ber. In particular, the Gutzwiller approximation as- 
sumes the ground state of the system has the form 

l*Mi.) = nfil {Y.n=ofn I"),), where the f^'^ is the 
amplitude of finding n particles on site i [2^. The co- 
efficients fn^ can be determined by variation to mini- 
mize the total energy {"^ mf\ Hbh |^'a/_f)- In contrast to 
the state of Eq. ([2]), |^'a/_f) is neither an fixed- number 
state nor a fixed-number state (i.e., its total number 
fluctuation AA^^ jg nonzero; AA^^ = ^/^2^ _ ^jy)2 and 
A'^ = X^i^-i)- Assuming a Gutzwiller state, by substi- 
tuting a\aj = (al)aj + (ai)al ^ into the Bose- 
Hubbard Hamiltonian and diagonalizing it, one obtains 
the mean-field excited states [26| . 

In the deep lattice regime, when [{n — 1)U + nZJ] < 
fi < [nil — (n 4- 1)ZJ], where Z is the coordination num- 
ber of the lattice, the Gutzwiller state is ||n)) -I- 0{J/U) 
with zero condensate fraction. In this range of the chem- 
ical potential, the system has zero compressibility (de- 
fined as the change in the expectation value of the total 
number of particles with respect to the chemical poten- 
tial or d{N)/dfj.). This state is identified as the Mott- 
insulator. When [nU-{n+l)ZJ] < ^ < [nU+{n+l)ZJ], 
the Gutzwiller state is l\i\ip)i + 0{J/U) with = 
fn \n) + fn+i \n+l). Because each single-site state is a 
linear combination of only two states, the system can be 
interpreted in the l anguag e of spin physics (the pseudo- 
spin treatment (l5l.l2ll. [27lf28jl. The single-site state can 
be rewritten as 



IV') = e'^/^ sin-\n)+e-'^'^ co&-\n + l)+0{J/U), (4) 



where the parameters 6 and (j) represent the spherical 
angles of the expectation value of the spin. In the 
pseudo-spin model, the Mott-insulator state is repre- 
sented by a state where the pseudo-spins are either all 
up or all down along the z direction. The condensate 
state is represented by a state where the pseudo-spins 
point in a direction determined by the chemical poten- 
tial and a 4> direction determined by spontaneous sym- 
metry breaking. We set (/) = for convenience and find 
that the relation between 9 and the chemical potential 
is cosO = {fj, — nU)/[{n + 1)ZJ]. The number fiuctua- 
tions on a single site, the density matrix elements, and 
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100 sites, long. The parameters of the optical lattice are 
A = 1.064 and Vo = 20i?/f and the frequency of the 
external harmonic trap is cj^ = 27r x 6.13 Hz. The system 
has an n = 2 Mott-insulating sphere, an n = 1 Mott- 
insulating shell, a condensate shell between the two Mott- 
insulators, and a condensate shell on the outer boundary. 
The thickness of the condensate interlayer is 2.66/im or 
5 sites long, approximately ten times thinner than the 
Mott-insulating regime. We use this setting to simulate 
the rf spectroscopy experiment in Sec. 111. 



FIG. 1: (Color online) Number density profile (solid line; scale 
on left axis of graph) and condensate fraction (dashed line; 
scale on right axis of graph) of an inhomogeneous lattice bo- 
son system. From inner to outer of the system, the phases 
(thickness in lattice spacing ds) are the Mott-insulator (55), 
condensate (5), the Mott-insulator (38), and condensate (2). 

the condensate fraction when < 6* < tt are found to be: 
An^ = cos2(6i/2)sin2(6i/2), 

{a\aj^i) = (n + l)An^, 

fc = {a\aj^i) / {a\a,) . (5) 

By setting the chemical potential to the value which 
renders the expectation value of the total number of par- 
ticles of the mean-field ground state equal to the total 
number of particles of the fixed-number ground state, we 
find the relation cos^(0/2) = Mi/M. In turn, this re- 
lationship equates the density matrix of the two states 
given by Eq. ^ and Eq. (O, thus showing that these 
two states concur in crucial aspects of the condensate 
and justifying the usage of the mean-field state for cal- 
culational purposes. In our subsequent calculations, we 
therefore use whichever representation of the condensate 
phase is most convenient. 

C. Inhomogeneous system 

For inhomogeneous lattices, local condensate order can 
be studied by replacing the chemical potential ^ with 
(/i — Vi), corresponding to the local density approxima- 
tion. Because of the spatial variance of Vi, the sys- 
tem is expected to display Mott-insulator phases and 
condensate phases that are spatially separated [1, 
In a three-dimensional (3D) deep lattice, for example, 
if the potential is a radially symmetric harmonic trap 
{Vi = mijjlr'l /2), the system will have a concentric 
shell structure with thick Mott-insulator shells alternat- 
ing with thin condensate shells. The density profile of 
this system resembles the side view of a wedding cake 
(Fig.©. 

In Fig. [U we plot the profile for an experimentally 
feasible system containing = 5.11 x 10^ *^Rb atoms 
with as — 5.32 nm in a sphere whose radius is 53.2 /im, or 



III. RF SPECTROSCOPY 

One effective probe of the structure of lattice bosons 
is the high-resolution microwave spectroscopy. In re- 
cent experiments, an external rf is applied to in- 
duce ®^Rb atoms to make hyperfine transitions from 
|F = l,mF = -1) state to |F = 2,tof = 1) state 
or to |F = 2, mi? = 0) state [l^l- In this section, we pro- 
pose variants of the experiment of Ref. (Tsj in which the 
Mott-insulating shells with different occupation numbers 
have been resolved by analyzing the density-dependent 
clock shifts of the system. The resolution hinges on the 
energy differences of interaction among two atoms in the 
same hyperfine state versus in different hyperfine states. 

In order to probe the condensate in the deep-lattice 
regime, we propose specific values of lattice parame- 
ters that would distinguish the spectrum of the Mott- 
insulator state and that of the condensate state in two 
ways. The first is that the former would have one reso- 
nant peak, but the latter would have two. The second 
is that the two peaks of the latter would be blue-shifted 
with respect to that of the former. The scales for the 
separation between the two peaks and the shift are of 
a order of a few Hz, and we anticipate that they will 
be within experimental reach in the near future, pro- 
viding an effective means to verify the existence of the 
condensed interlayers. By considering the effect of Gold- 
stone modes and finite-temperature, we show that this 
signal is a signature of condensate order, and would be 
absent if the interlayers were a normal boson fluid of the 
same density. 

Below, we begin with a discussion of our rf setup ap- 
plied to a homogeneous system and then use these re- 
sults to derive the spectrum of the inhomogeneous sys- 
tem in the local-density approximation. We use nu- 
merical simulations where appropriate and also invoke 
Fermi's golden rule to analyze the rf spectra for a range 
of phase space. To summarize our findings: at zero- 
temperature, the Mott-insulator and condensate states 
can be clearly distinguished via their single-peak versus 
double-peak structures. For the set of proposed parame- 
ters, rf transitions into states containing Goldstone mode 
excitations (which are the low-energy excitations asso- 
ciate with the condensate) do not obscure the double- 
peak structure at zero temperature. However, at fi- 
nite temperature, thermal excitations affect the two-peak 
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Chemical Potential (yu) 

FIG. 2: (Color online) Energies of the lowest four single-site 
states of the 2-species Bose-Hubbard Hamiltonian. The corre- 
sponding states, which are defined in Eqs. (ffl Ulf) . are marked 
for each curve. In the mid-region (marked red), the ground 
state of the system is a condensate, while on the sides (marked 
blue), it is a Mott-insulator.In each region, the ground state 
can make transitions to those higher energy states denoted 
with solid lines (with corresponding energy gaps marked by 
the black arrows), but is forbidden from transitions to those 
states marked with dashed lines. 

structure. At very low temperatures (ksT <C ZJ), Gold- 
stone modes lead to a temperature-dependent broadening 
of the peaks in the condensate signal. At higher temper- 
atures {ksT ^ ZJ <C U), thermal fluctuations destroy 
the condensate order to yield a normal fluid; ksTc is of 
order ZJ [ll,[33|. In this temperature regime, Goldstone 
modes completely obliterate the peaks, and extra struc- 
ture develops at other eharacteristic frequencies corre- 
sponding to new allowed hyperfine transitions. At tem- 
peratures ksT > U, the system displays large density 
fluctuations, even in the regions that are Mott-insulating 
at zero temperature, and the system loses signs of its 
quantum phases. Our analysis is thus limited to tem- 
peratures much less than the interparticle interaction, 
[/ ~ fcs X 10 uK wherein, as a function of temperature, 
the resolution of the peaks in the rf spectrum tracks the 
condensate order and its ultimate destruction at Tc. 

A. Zero-temperature rf spectrum: Decoupled-sites 

For a uniform system of two-state (|a) and \b)) bosons, 
the Bose-Hubbard Hamiltonian contains two tunneling 
strengths (Ja and Jb), three interaction strengths (C/q, 
Ub, and Uab), the chemical potential, and the single- 
particle energy difference between a and b particles (wba)- 
In the following, we consider a particular case: Jb = 0, 
Ub^ Ua, Uab = Ua, and ojba ^ Ua- This setting has two 
advantages. (1) Because ujba > 0, the ground state of the 
system has all particles in the a hyperfine state but no 
particles in the b state. (2) Because Ub ^ Uab — Ua, the 
energy gap between having two b particles on a site and 
having one b particle is much larger than that between 



FIG. 3: (Color online) Absolute value of transition amplitude 
Alt and A2t (solid and dashed lines, respectively; axis on 
left-hand side of graph) and condensate fraction (dotted line; 
axis on right-hand side of graph) as a function of the value 
of the cosine of the pseudo-spin angle (corresponding to the 
density of bosons) . The system is a uniform condensate whose 
number density satisfies 1 < (n) < 2. 

having one b particle on a site and having no b parti- 
cles. Therefore when the frequency of the rf field is of 
the order of the gap in the latter case, we can safely limit 
consideration to transitions to states with one b particle. 

For the situation described above, the Bose-Hubbard 
Hamiltonian, including the upper hyperfine state, be- 
comes 

-J ^ {aluj + d^CLi) + '^[^{flai + nbi){nar + flbi - 1) 
<ij> i 

-Kf^ai + flbi) + ^ ^ flbi{nu - 1) + i^baflbi], (6) 

where we drop the subscript a of J a and C/„. In the 
deep lattice regime, we use the mean-field approximation 
to analyze the rf spectrum of the system. The single- 
site Mott-insulator state with (n) = n and the single-site 
condensate state with n < (n) < n+1 can be represented 
by 

IV'mW) = \n)<E)\Q), 

9 9 
\il>c{n)) = sin- |n) (K) |0) +COS- |n-t- 1) ® |0) , (7) 

where \na) ® \nb) represents a single-site state with 
particles of hyperfine state a and Ub particles of hy- 
perfine state b. As in Eq.® of See. IIIBl we have 
COS0 = (/i — nU)/[{n + ^)ZJ] determining the angle 
between the pseudo-spin and the z axis. Compared to 
Eq.dS]), the symmetry-breaking phase has been set to 
zero for the decoupled-site analysis; the subsequent Gold- 
stone mode analysis implicitly assumes that this phase 
has gradual variations from site to site. 

The Hamiltonian describing the interaction between 
the bosons and an applied rf field is Q 

Hrf = ^ l{a^ibi + i>\oii) cosujt, (8) 

i 

where 7 is proportional to the amplitude of the rf field 
and uj is its frequency; this Hamiltonian can be de- 



6 



rived from the second-quantized form of the interaction 

We analyze the allowed transitions in different ranges 
of the chemical potential. We evaluate the transition 
plitudes At = (^IX]i7(^l^j + \^) between allowed <-> 
initial and final states |/) and \F), respectively, and the 
energy gaps AE, given by the difference in the expec- 
tation values of the Bose-Hubbard Hamiltonian of Eq. 
([S]) in the corresponding excited state and ground state. 
We plot the energy levels in Fig. [2l which shows how the 
system at zero temperature undergoes quantum phase 
transitions through the n Mott-insulator state [\tjj„i{n))], 
the condensed state [|-0c('t-))]! and the n -I- 1 Mott state 
[\ipm{n + 1))] as the chemical potential is increased. The 
figure also shows possible excited states to which the 
ground state can make a transition in the presence of the 
rf field and shows the energy gap between them. Details 
of the transitions for each range of the chemical potential 
in Fig. [5] are discussed below. 

When [{n - 1)U + nZJ] < fi< [nU - {n + 1)ZJ], the 
ground state is the n Mott-insulator state [|V'm(fi))]. We 
have only one possible excited state [denoted by |em(n))] 
to which the ground state can make a transition. The 
excited state, the transition amplitude At, and the energy 
gap AE are correspondingly 
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FIG. 4: Fraction of bosons transferred to the higher hyperfine 
state for a uniform condensate with 1 < (n) < 2 at t = 100, 
400, 1000, and 2000 ms. The horizontal axis is the detuning of 
the rf field with respect to the resonant frequency of the Mott- 
insulating state, while the vertical axis is cos 6. The color scale 
is bounded by (black) and 1 (white) with contours separated 
by 0.1. 



Because AE is independent of n, all the Mott-insulator 
states with different n have the same energy gap between 
the ground state and the excited state. 

When [nU - (n + l)ZJ] < < [nU + {n + 1)ZJ], the 
ground state is the condensate state \ipc{n)). We find 
two orthogonal excited states [denoted by |eci(n)) and 
|ec2('T-))] with nonzero transition amplitudes. The first 
state, its corresponding transition amplitude Ati, and 
energy gap AEi are 



AE2 are 

6 9 

\ec2{n)) = cos y |?i - 1) (g) |1) - siny |n) (g) |1) , 

^2t = 7(V^sin^cosy - Vn + Icos^siuy), 

AE2 = ujba + ^'^^^^^ [V^TTlism^ 6 + 2) 

-I- sin 6* sin 6*1 + + 1 cos 6' cos 6'i)](ll) 



|eci(n)) = siny |n- 1) (g) |1) +cosy |n) (g) |1) , 
A\t = 7(V'^sm-smy + Icos-cosy), 

Ai?i = u;ba + ^'^^^^ [V^(sin^ B + 2) 

- sin 6* sin 6*1 + 2V?i + 1 cos 6* cos 6*1)] (10) 

where [ed) is taken to be the equilibrium state (the 
lowest energy state) of the Hamiltonian with the con- 
straint Ub — 1 (i.e., in the nb = 1 block). We calculate 
the parameter 61 by minimizing the energy of the sys- 
tem with site i in an excited state [|eci(n))] and all the 
other sites still in the ground state [|i/'c(?t-))]. As a re- 
sult, we obtain a relation between 9i and 9, which is 
tan6'i = ^yn/{n+ 1) tanfl. The second excited state, 
corresponding transition amplitude At2, and energy gap 



When the chemical potential in the condensed phase 
reaches its upper or lower bound (6* = or 6* = tt), Eq. 
([T0|) becomes Eq. ([9]) with n or n+1 Mott-insulator states 
correspondingly, but the transition amplitude of Eq. (jlip 
vanishes (Fig. [3]). However, the n = condensate state 
is an exception since it can only make a transition to one 
excited state |0) g) |1) with a transition amplitude 7 cos | 
and an energy gap {uba + J^sin'^ |). 

Thus, while the Mott state can only make a transition 
to one excAted state, the condensate state can make tran- 
sitions to two excited states. 

As a consequence of these allowed transitions, the 
Mott-insulator state ought to have a single peak in the 
associated rf spectrum while the condensate ought to 
have two. This can be seen to first order (in number of 
particles transferred) by employing Fermi's golden rule 
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Average fraction transferred 




FIG. 5: Average fraction of bosons transferred between t = 
1000 and 2000 ms. Contour interval is 0.1. The first peak 
appears in the range of to 7 J (6.29 Hz), and the second 
peak appears at about 11. 5J (10.3 Hz). When cosO = ±1, we 
recover the Mott signature. 



(FGR); where the transition rate is given by [3l| 



^(^) ='^Y.^Pi- 1^*1''^ {lo~Ep + Ej) , (12) 

FJ 

where |/) {\F)) is the initial (final) state with energy 
Ei {Ef) and probability for occurrence pi (pp)- Hence, 
FGR would predict a single delta function peak for the 
Mott-insulator state rf spectrum and two delta function 
peaks for the condensate. 

We now confirm that the presence of the condensate 
does yield a two-peak rf spectrum by numerically time 
evolving the equations of motion in the presence of the 
rf field. Taking into account both experimental feasi- 
bility and optimal conditions for observing the rf signa- 
tures of the condensate, we consider the setting of Fig. [T] 
This setting corresponds to the values U ^ 2tt x 306 
Hz, J ^ 2n X 0.898 Hz (then ZJ/U = 1.76 x lO^^ is 
small compared to 1), and p = 1.5U. The intensity of 
the rf field is set such that 7 = 0.118ZJ = 2n x 0.637 
Hz, which is of the same order as that in the experi- 
ment of Ref. If the rf field is turned on at t = 0, 
within the rotating-wave approximation we find that the 
fraction of transferred particles at the frequency of the 
second peak of the condensate's spectrum begins to rise 
at t = 0.225/7 (400 ms), and that at any given time, 
the spectrum shows oscillatory behavior with respect to 
the detuning (Fig. |4]). In order to see a distinct second 
peak, we calculate the average fraction transferred be- 
tween t = 0.637/7 (1000 ms) and t = I.27/7 (2000 ms). 
Figure [5] shows the average fraction transferred for a uni- 
form condensate. We can see two peaks whose positions 
agree with the energy gap in Eqs. pO]) and (fTTj) . Figure 
|6] shows the average fraction transferred for the entire in- 
homogeneous system. Because the fraction transferred in 




-5 5 
Av(J/h) 




FIG. 6; (Color online) Average fraction of bosons transferred 
between t = 1000 ms and 2000 ms for the inhomogeneous 
system defined in Fig. [1] Up: the system with both Mott- 
insulator and condensate regions (solid) compared with that 
having only Mott-insulator regions (dotted). Down: spec- 
trum of the condensate shells after eliminating the contribu- 
tion from the Mott -insulator shells. The results show that 
two peaks with positive shifts are signatures of condensate 
order in a deep optical lattice. 



the Mott-insulator is an even function of the detuning, we 
subtract the negative detuning part from the correspond- 
ing positive detuning part to eliminate the contribution 
of the Mott-insulator. The result shows that there are 
two peaks of magnitude 1% (5 x 10^ atoms transferred), 
and their positions agree with our theoretical analysis. 

All calculations above rely on a mean-field approxima- 
tion. As a minor check that number conservation and en- 
tanglement do not alter the proposed signatures, we per- 
formed the numerical toy simulation os time-evolving the 
entangled, fixed-number condensate state of six bosons in 
four sites in the presence of an rf field (Appendix C). Even 
such small system sizes has been shown to exhibit salient 
features of the superfluid-Mott-insulator phase transition 
in macroscopic systems [32|. The results show the same 
signatures as those of the mean- field approach: two peaks 
with positive shifts of order J, thus providing some con- 
fidence that our result is not an artifact of the mean-field 
approximation. 

To ensure that no spurious effects emerge from truncat- 
ing the Hilbert space, we have analyzed the effect of the 
leading subdominant states \n — 1)®|0) and \n + 2)^\0). 
Within the resultant four-state truncated space, we ob- 
tain four possible final states. Two of them have the 
same energy shifts (of order J) and transition ampli- 
tudes (of order 7) as those obtained in the two-state trun- 
cated space, plus small corrections of order J/U; these 
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corrections contribute a fraction {J/UY to the number 
transferred. The other two final states have energy shifts 
of order U J) and transition amphtudes of order 
"fJ/U 7), which also contributes a fraction {J/UY' 
to the number transferred. The minimal change coming 
from taking into account these additional states justifies 
our usage of the two-state truncation in the deep lattice 
regime. 



B. Goldstone modes and finite-temperature effects 

At finite temperature, the thermally-occupied excited 
states of the system will affect the rf spectrum. When 
the temperature is much smaller than the interaction en- 
ergy {ksT <Si U), the Mott-insulator state has negligible 
thermal effects because it has an excitation gap of order 
U [1, m, [l!]. However, we should be concerned that the 
two-peak signature of the condensate may be destroyed 
at low temperatures by Goldstone modes associated with 
the continuous symmetry-breaking in the phase (j) of the 
condensate ground-state wave function. Each Goldstone 
mode or boson corresponds to a long-wavelength dis- 
tortion of the density and phase between neighboring 
sites and the modes form the gapless low-energy exci- 
tation spectrum for the condensate. The structure of 
these modes has been derived from the pseudo-spin for- 
mulation and analyzed for homogeneous and inhomoge- 
neous systems in Ref. [H]. In order to be self-contained, 
we present the Goldstone-mode description here in Ap- 
pendix B for the homogeneous generalized case of the 
two-species system of interest in this paper. 

At zero temperature, the initial state consists of the 
ground state of the system purely comprised of a parti- 
cles, and thus no Goldstone modes are excited. Under an 
rf field however, a small number of particles change their 
internal state to b. If this change is not uniform in space, 
it will be accompanied by Goldstone modes. Assuming 
that Ng is the total number of Goldstone bosons excited 
by the rf field, we can estimate the final-state energy 
of a single site, which differs from that of the Nq = 
case, as {Ng/M)ZJ, where M is the total number of 
sites and ZJ is the energy scale of the Goldstone-mode 
spectrum. If the rf field is weak enough that the average 
number of excited Goldstone bosons per site {No /M) is 
much smaller than one, the change in the single-site en- 
ergy gap is much smaller than ZJ, which is about the 
distance between the two peaks obtained in the Sec. II. 
In typical experimental settings, this indeed is the case 
given that about 10% of the particles make transitions to 
the b state. Therefore, exciting Goldstone modes at zero 
temperature ought not obscure the two-peak signature; 
in Appendix B, we explicitly show that these modes pro- 
vide a small background which still leaves the two-peak 
signature robust. 

At finite temperature, in three dimensions, the con- 
densate density is reduced by the existence of thermally 
excited Goldstone modes. As argued in Appendix B, for 



low temperatures kT ^ ZJ, we expect a slight reduc- 
tion in the height of the two peaks and a slightly larger 
contribution to the background. Nevertheless, for three- 
dimensional systems, we expect the two-peak condensate 
signature to persist at these temperatures. 

Close to the critical temperature kTc « ZJ (see Ref. 
0j d^l), we expect the average thermal expectation 
value '^j{o-j)/N of the order parameter to vanish, cor- 
responding to destruction of long-range order in the sys- 
tem. Within the decouplcd-site mean-field approxima- 
tion, {tti) fluctuates in magnitude and phase from site 
to site. The kinetic part of the mean-field energy in the 
Bose-Hubbard Hamiltonian can thus be thought of as 
having a spread of order (n + 1)ZJ. Given that the two 
low-temperature peaks in the rf spectrum are also sepa- 
rated by about (n+ 1)ZJ, we predict that the two peaks 
merge into a continuum around the critical temperature. 
In two-dimensions, even the smallest temperature suffices 
for the Goldstone modes to destroy true long-range order 
and we believe that this will be reflected in the smearing 
out of the two-peak structure for low-dimensional sys- 
tems. 

Finally, for the inhomogeneous shell situation, the 
Goldstone modes become quantized due to confine- 
ment [1^. Reference [HI provides a discussion of the 
mode structure and dimensional crossover as a function 
of temperature. The essential issue is whether the quan- 
tized modes along each direction which are accessible at 
temperature T are numerous enough to form an effec- 
tive continuum. Modifying the arguments above to in- 
clude these inhomogeneous effects, we deduce that the 
two peak structure would remain at low temperatures 
for a condensate shell with a thickness of several lattice 
sites. 

In summary, our arguments indicate that the presence 
of condensate order can be detected in a two-peak rf spec- 
trum for appropriate rf parameter settings in contrast to 
the one-peak structure of the Mott-insulator phase. Our 
results also suggest that when the condensate becomes a 
normal fluid at T = Tc, the two-peak structure is washed 
out. Our arguments can be made more rigorous by way 
of numerical simulations such as those of Sec. IIIA that 
would include the Goldstone modes and finite tempera- 
ture effects; such a treatment is beyond the scope of this 
paper. 



IV. MATTER WAVE INTERFERENCE 

Experimental observation of matter-wave interference 
peaks in absorption images have provided striking evi- 
dence for superfluid states in shallow optical lattices Q. 
In addition, the concentric Mott-insulating-shell system 
has been observed to preserve finite visibility in the in- 
terference pattern [H, [1^ . Characteristics of the inter- 
ference patterns have been analytically and numerically 
studied under several specific conditions for normal fluid, 
superfluid and Mott-insulating phases in the lattice bo- 
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son system [37|, |38|, |39|, |40|. Here, we consider bosons 
in the deep lattice regime and ask what signatures of 
matter- wave interference can distinguish the presence of 
a condensate shell between two Mott-insulating shells. 
We work at zero temperature, which should be valid for 
temperatures kT <C J, where one would expect the con- 
trast between these states to be strongest. We begin 
by considering a homogeneous system, where the many- 
body wave functions of the Mott-insulating and deep- 
lattice condensate phases can be exactly time-evolved to 
obtain the density profile measured by absorption imag- 
ing. For simplicit y, w e ignore the effect of interaction dur- 
ing time of flight (37[ , which is expected to quantitatively 
influence the intensity and width of peaks but preserve 
the qualitative signatures of the interference pattern [s^ . 
By calculating the time-evolution of a system of bosons 
on a ring lattice, we show that condensate states will have 
sharp interference maxima, in contrast to Mott-insulator 
states, after free expansion of the system. To illustrate 
how these features would be displayed in an inhomoge- 
neous system, we also the time evolution of a system 
of concentric ring lattices that have a Mott-condensate- 
Mott structure and point to signatures of the condensate 
in the time-evolved profiles. 



A. Density profile upon expansion for 
homogeneous systems 

For the Mott-insulator state | , the many-body wave 
function is a product of A'^ single-particle wave functions: 

«'(ri,...,rjv) = (ri,...,rjv | \n)) 

N 

sym{r} i— 1 
{s} i=l 

where — s^*^) is the single particle wave function 
of particle i localized on a site at position s^*^ and A 
is a normalization constant. To satisfy Bose-Einstein 
statistics, we need to symmetrize all N degrees of 
freedom; this symmetrization (denoted by X]sym{r}) 
is equivalent to the sum over all configurations of 
{s*^*'} (denoted by X]{s})j which are permutations of 
{ si, - y ,si , S2, ,S2 , . . . sa/,.^.,sa/ }, where Sj is the po- 

n n n 

sition of site j. 

Near the center of a lattice site, the lattice poten- 
tial can be approximated as harmonic with frequency 
u) ~ \/4VoEfj/h, where Vb is the depth of the optical 
lattice and Er is the recoil energy. The single-particle 
wave function can thus be approximated as a harmonic 
oscillator ground-state wave function well-localized on a 
site. If we ignore interactions (valid at low density), af- 
ter turning off the lattice potential and the trap potential 



all single-particle wave functions are time-evolved by the 
free particle Hamiltonian p^/2m. The wave function at 
later times is then given by 



-sW) 



2/2 \ "^/^ exp[-(r, - sW)V(/2 + 2iht/m)] 
(P + 2iht/my/^ ' 



(14) 



where d is the dimensionality of the system, and I is 
the characteristic length of the single-particle wave func- 
tion [13, EH . The density profile as a function of position 
and time is obtained by calculating the diagonal elements 
of the single-particle density matrix [23| ; 

p{r,t) = iVy"^'j(r, r2, . . . ,rAr)'I't(r, r2, . . . ,rjv)(ir2 . . .drAT 



N\AfY.Y.^Ur~s('^)Mr-s'^'^) 

{s} {s'} 



N 



exp[-^(s« -s'^^^)V(2/2^ 



(15) 



For the condensate state \'i>) of Eq. ([2]), the density 
profile is given by 

Pe(r,t) = N\AfJ2Y.^W}^M 

M W} 

X { E E ^:(r-sW)^t(r-s'(^)) 



{s;M}{s';{r,'}} 
N 

exp[-^(sW-s'«)V(2/^)]}, 



(16) 



where {»;} denotes the same sum in Eq. and 
S {s:{))}} denotes the sum over all {s*-'^} with a specific 
{77}; {s*^*^} are permutations of the set including n + 1 
each of the sites 771 , ... , ?7j\/i and n each Sj of the 
other sites. If Mi = or Mi = M, {?/} has only one 
configuration, and hence Eq. becomes Eq. 



B. Example: Expansion of bosons in a ring lattice 

As a concrete case that shows how sharp interference 
peaks emerge upon release and expansion of the conden- 
sate, we consider the illustrative example of bosons in a 
one-dimensional lattice of sites located on a ring of radius 
s in a two dimensional plane. Because of the symmetry of 
the ring, constructive interference is expected at its cen- 
ter. We therefore calculate the center density contributed 
by one particle, defined as p{t) = p{0,t)/N. The result 



~u. ,'2l'V^' exp[-2.2p/(/4 + Ah^tym^)] 
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where d ^ 2 for this system and F is equal to 

N 

{s} {s'} ^=2 

for the Mott-insulator state and 
{')} W} 

N 

x{ E E exp[-^(s«-s'«)V(2/^)]} 

{s;M}{s';{r,'}} »=2 

for the condensate. If s 3> ^ p{t) has a maximum at 

t„, mls/Vdh (18) 

and a half width At = 2^/'^rals/h. While tm indicates 
the time it takes matter waves to disperse from the ring 
to the center, At measures the duration of constructive 
interference processes. Equation (|18p is of the same form 
as a continuum superfluid confined to a shell geometry 
with the characteristic length I given by the thickness of 
the shell [l^]. The terms with {rj'} ~ {77} in the sum for 
F are the same for the Mott-insulator and condensate 
states. However, F for the condensate state has many 
additional terms with {77'} ^ {77} that contribute to the 
center density which are not present in Ff,[. Therefore, 
we see that the center density of the condensate state 
has a much sharper peak than that of the Mott-insulator 
state. At finite temperature, the interference pattern is 
the thermal average over the densities of all possible pure 
states. The coefficients C{^} of higher energy states are 
not necessarily real and positive and thus would decrease 
the value of F given that not all terms give positive con- 
tributions. Therefore, we expect that the interference 
pattern of a condensate becomes more blurred with in- 
creasing temperature. 

Although here we have only analytically solved for the 
density at the center of the ring, constructive interference 
patterns will occur throughout space after the system 
is released from its ring trap. The differences in these 
spatial patterns between initial Mott-insulator and initial 
condensed states gives further evidence for condensate 
order. As a toy example that exhibits these patterns, 
we numerically simulated the expansion processes of a 
four-site Mott-insulator and a four-site condensate in the 
ring geometry (Appendix C). Although the initial density 
profiles of the two cases look the same, their interference 
patterns behave quite differently during the expansion 
process. 

C. Inhomogeneous systems 

The above discussions show how homogeneous systems 
of condensate and Mott-insulator display different pat- 
terns upon release and expansion and that interference is 




0.0 0.5 1.0 1.5 2.0 

t(lO-Vc) 

FIG. 7: (Color online) Density at the center point, as a func- 
tion of time, during expansion of the system. The inset shows 
the geometry and the site positions of the model. The solid 
curve corresponds to the system that has a center site and an 
outer ring in Mott-insulating states and the inner ring in the 
condensate state. The dashed curve corresponds to the con- 
trol system with the inner ring in the Mott-insulating state. 
The vertical axis is the relative density with respect to the 
initial value. The peak and the relatively large value at long 
times are evidence of the condensate part of the system. 



a signature of condensate order in the system. To illus- 
trate how these signatures can be used to probe conden- 
sate order realistic experimental situations, we consider a 
simplified inhomogeneous system of ®'^Rb atoms display- 
ing the nested, or "wedding-cake," structure of bosons in 
optical lattices in the presence of a harmonic trap. Build- 
ing on the results of the previous section for bosons on 
a ring lattice, we model a two-dimensional system made 
up of a center point and two concentric rings, where the 
central site is in the n = 2 Mott-insulating state, the 
four sites in the inner ring are in the (n) = 1.5 conden- 
sate state, and six sites in the outer ring are in the n = 1 
Mott-insulating state. Site spacings of the two rings are 
both ds = 0.532 fim. For the lattice depth Vo = 20Er, 
the characteristic length of the initial single-particle wave 
function is / — 0.283ds. If the lattice potential is turned 
off at t = 0, the density at the center point as a func- 
tion of time shows a peak at t = 4.2 x 10~^ s (Fig. [7]). 
Compared to a control system where the inner ring has 
Mott-insulator states of the same density (2 sites in n = 2 
and 2 sites in ti = 1), the peak and the larger long-time 
residual density in the center indicate large constructive 
interference due to the condensate part of the system. 
Considering expansion of a 3D inhomogeneous system of 
Fig. [1] we expect to observe a density peak in the center 
of the sphere at < = 3.67 ms, which is about 1/4 the du- 
ration of the time of flight in the experiment of Ref. (3] . 
We expect that such a time-evolved density profile should 
also be able to discern the presence of condensate inter- 
layers for larger systems with realistic numbers of bosons. 

In summary, we obtain interference signatures of a thin 
condensate inter-layer between two Mott-insulator layers. 
For a toy system exhibiting concentric co-existent phases, 
we show that the condensed phase contributes a sharp 
peak to the time-evolved central particle density. As per- 
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formed for single phases, simulations of realistic inhomo- 
geneous large-sized three-dimensional systems within our 
truncated wave function basis and extraction of measur- 
able quantities, such as the visibility, are in order; we 
believe that our results here are indicative of the signa- 
tures of the condensate layer that would be obtained in 
actual systems. 



V. SUMMARY 

Systems of atoms in optical lattices provide a 
seemingly-ideal system for exploring the quantum phases 
of bosons. Generically, experimentally-studied systems 
are not of uniform density, due to the confining trap. This 
leads to the coexistence of Mott-insulating and conden- 
sate phases within a single system at fixed lattice depth 
and complicates the study of these phases. Here, we have 
considered bosons in deep optical lattices, where conden- 
sate phases exist between Mott-insulator phases, but arc 
close to their quantum transition to the Mott-insulator 
phase. We have justified the use of a mean-field approx- 
imation by constructing the entangled, fixed-number, 
many-body wave function and showing that many key as- 
pects of this phase (in particular, the condensate fraction 
and the on-site number fluctuations) are identical in the 
two descriptions. Toward developing experimental tests 
of the condensate inter layers in inhomogencous systems, 
we have calculated the condensate signal using two com- 
mon probes of cold-atom systems: rf spectroscopy and 
release-and-expansion. We have shown specific features 
of the signal in these probes that distinguish the conden- 
sate regions from the Mott phases. We have shown that 
the rf signal specifically probes condensate order in that 
these characteristic features disappear when the conden- 
sate is replaced by a normal fluid as a result of thermal 
fluctuations. We hope that our work will enable closer 
exploration of the structure of coexisting Mott-insulator- 
condensate phases in realistic experimental situations. 
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APPENDIX A: DENSITY MATRIX OF THE 
FIXED-NUMBER CONDENSATE STATE 

In this appendix, the elements of the density matrix of 
Eq. @ arc calculated. For diagonal terms. 



Ml 

n 



Ml 

n 



" E iQr,}!' + (« + !) E I^mI' 
" + E I'^wl'' 



where {(x\ai)i 



i {n\alai\n)^. From symmetry, 

X]ie{r;} l^fi/lT independent of i but depends only on 
the number of terms in the sum. Considering the nor- 
malization condition, we find 

El 1 2 # of terms in the sum 
\'-'{ri}\ = 



ie{r,} 



^ of terms in the sum ^ 

M -l\ 
Ml -I ) _ Ml 

~ ^ ^ m" 



{')} 



{a\ai) = n + 



M 
Ml 

Ml 

m" 



(Al) 



To evaluate the number fluctuation of a single site, we 
need to calculate (a|aiajai). Similarly, 

{alaialai) = + {2n + 1) ^ |C{r,}|^ 







ie{v} 


+ {2n 




Ml 

'"m 


(alcLiala,, 






M ^ 


Ml 
M ' 





(A2) 



For the off-diagonal terms 

(n + l) 



{l)2---')A/i } 

where {772 .. . r]M^ } is a set of distinct integers chosen 
from {1,2,...,M} excluding i and j and J2 {v2---vmi} 
denotes the sum over all combinatorial configurations. 
Because C^^j^ is determined by the energy cost of the 
corresponding state, the difference between jj^.. .^^^^ j 



12 



and C{j,T;2...j)Afi } estimated to be of order 1/Mi, and 
hence vanishes in the thermodynamic Hmit. By consid- 
ering symmetry and the normaHzation condition, we find 



a a ,- 



(77, + 1) ^ |C{i^,,2...^j^,J 
{ri2---riMi } 

# of terms in the sum J2 {v2---riM, } 

1 ) ■ Z 

^ of terms in the sum ^ {^j. 
M -2 

(77+1)^^ ^ = (?7+ 1)A772. (A3) 

Ml 



In the position representation, the density matrix is 
of the form with all diagonal elements x = n + Mi/M 
and all off-diagonal ones y = {n + 1)A77'^ and can be 
diagonalized by Fourier transformation to the momentum 
representation: 



{x-y)I + y^\i) (j| 

{x-y)Y,\k) {k\+yM\k^Q) (fc = 0| 



(A'/ - l)y] \k = 0) (fc = 0| + (x - y) ^ |fc) 

fe#0 



(A4) 



Because M is much larger than 1, the fc = state has 
much larger occupation than all the other states, which 
are uniformly occupied in a small fraction. The conden- 
sate fraction of the system is given by [x + {M — l)y]/N. 



APPENDIX B: PSEUDO-SPIN MODEL AND 
GOLDSTONE MODES 



In the rf transition process, the system with M sites 
can be well described in the mean-field approach as a 
product of M single-site states. Up to leading order in 
J/U, each single-site state is represented within a trun- 
cated basis of four vectors, |77)(3)|0), |n-|- 1) (g) |0), |77)(3)|1), 
and |77 + 1) (g) 1 1) , where wc use the notation \na + nin rih) . 
As each number has only two possible values in the trun- 
cated bases, we can map these states onto the states of 
two spin-1/2 spins as U) ® Ub), | t> ® life), U) ® | Tb), 
and I t)<8'| Tb) respectively. Consistent with this mapping 
of the basis, any operator in the Bose-Hubbard Hamilto- 
nian of Eq. ^ can be replaced by its corresponding spin 
operator as long as they have the same matrix represen- 



tation in the truncated space. Therefore wc have 

1 



77fa 



2 



1/2 -5f 5" 



1/2 + 5^5, 



b ' 



(Bl) 



where S = (7/2, and & is the vector of Pauli spin matri- 
ces. With these substitutions, we obtain the pseudo-spin 
Hamiltonian for the two-species Bosc-Hubbard Hamilto- 
nian of Eq. ([H]), up to a constant: 



Hps = -jY. V " + - V " + " 



S^SJ + SfS]) + ^ [{Un - t,)S: + . (B2) 



This Hamiltonian describes a system with two spins (S 
and Sb) on each site, each of which responds to an in- 
dependent external "magnetic field" in the z-direction. 
The interaction between transverse components of ad- 
jacent S spins has a strength related to the z- compo- 
nents of the b spins on those sites. When 3^,^ = —1/2, 
which corresponds to having no particles in the b hy- 
perfine state. Hps becomes the pseudospin approxima- 
tion to the single-species Bose-Hubbard model [lal- The 
ground state of Hps is the same as Eq. (O , when written 
in the corresponding spin basis. 

The equivalence of the matrix representations of the 
Bosc-Hubbard Hamiltonian and the pseudo spin Hamil- 
tonian imply that they also have the same excitations in 
the truncated space. In fact, the Goldstone excitations of 
the pseudo spin model, the spin waves, are equivalent to 
those of the Bose-Hubbard model, which are small den- 
sity and phase distortion over a long length scale. To 
obtain the dispersion relations of the spin waves, we in- 
troduce two species of Holstein-Primakoff bosons (A and 
B) to represent small fluctuations of the corresponding 
spins around their equilibrium value at zero temperature. 
Provided the fluctuations are sufficiently small such that 
all terms higher than quadratic can be neglected, the 
spin operators can be represented as S = {^^{A + A^) + 
sin 6»(i-iti)^|(i+it)^_ sine (i_^^t)+cos6i(i-iti)) 

and Sfc = (-i(B + B-f), ^{-B + fit), _ i + b^b). Af- 
ter substituting these into the pseudo-spin Hamiltonian 
of Eq. (|B2p and keeping terms up to quadratic order 
in the Holstein-Primakoff bosons, wc Fourier transform 
the resulting Hamiltonian into momentum space and per- 
form a Bogoliubov transformation on the A bosons into a 
bosons {Ap = UpCtp + WpciLp) to diagonalize the A terms. 
Thus we arrive at the diagonalized Hamiltonian: 



H 



PS 



■k/I 

= E 

p——'K /I 



BlBp, 



(B3) 
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N=4, t=0. 



N=2, t=0. 



0.00 0.05 0.10 0.15 0.20 0.25 0.30 
Detuning(U) 

FIG. 8: Average fraction transferred of a 6-boson 4-site en- 
tangled state driven by the external rf field. The parameters 
are Ua ~ Uab = 1, ZJ = 0.05, ^ = 1.5, and Ub ~ i^ta = 5. 
We can see the spectrum has 2 peaks with positive shifts of 
order J. 



where p is the lattice momentum and I the lattice spac- 
ing. The B terms correspond to the creation of a i? 
boson with momentum p, while the ap operator cre- 
ates a Goldstone excitation in the condensate with mo- 
mentum p. The excitation energy of the a excitations 



IS e„ 



zj{n + i)vTr 



Ip COS- 



where Jr, 



^ d ^cospil for d dimension. 

1=1 



For small p {pi <ti I) 



Ep = J{n -I- 1) Y Z sin 6{pl)^ + cos^ 6{pl)'^, which means 
that in the deep condensate regime {sin 9 ^ pi) the 
low-energy excitations are wave-like (energy oc p), while 
near the Mott-insulator boundary {sin 9 pi) they are 
particle- like (energy ex p^). The excitation energy of of 
the B particles is ujb = ^ba + (1/8)2' J sin 0, which is in- 
dependent of p. This excitation is gapped (reflecting the 
fact that the h bosons are in their Mott-insulator phase) 
and therefore should not be present in the initial state 
in rf transitions for temperatures kT <^ hcuba, which we 
assume throughout. 

Turning to the role of the Goldstone modes in rf tran- 
sitions, the form of the rf perturbation (given in the po- 
sition basis of the physical particles as ^^70^6^ -I- h.c.) 
is given in the Goldstone basis by 



t=2. 



t=2. 




t=4. 



t=4. 





t=20. 



t=20. 




FIG. 9: Matter wave interference patterns (per particle con- 
tribution). The left column is the 4-particle 4-site Mott- 
insulator, and the right one is the 2-particle 4-site conden- 
sate. The four rows from top to bottom are ai t — 0, 2, 4, 
and 20 time units(m/^/2/i). All parts of the figure have the 
same spatial scale. The number in the last row is the relative 
value of the density. 



Hrf - - 



j 7 sm ( 



E 



H.c. 



{vp + Up)Bl{ap + a^_p) 
(B4) 



where N is the total number of particles and n = n + 
cos^(0/2) is the number density at zero temperature. The 
first term creates a B boson uniformly in space without 
any distortions in density or phase (which arc represented 
by a bosons). The second term creates a B boson with 
momentum p accompanied by an a boson with momen- 
tum — p or annihilation of an a boson with momentum 
P- 

For the 3D system at zero temperature, using the 
above Goldstone representation of the rf perturbation in 



Fermi's golden rule [(Eq. [T^ ]. we find the transition rate 
due to the first term of Eq. (|B4[) is l'^^\uj) oc 5{uj — ujb)- 
This delta peak corresponds to the first peak we obtained 
in Sec. IIIA, although slightly shifted given the slightly 
different mean-field energy estimates. The transition rate 
due to the second term is /'^^^(w) oc sin^ 9\lo — lob\^ which 
is so small compared to the delta function near lob that 
the effect of J*^^) on the contrast of J*^^^ can be ignored. 
We expect similar argument to hold for the second peak 
obtained in Sec. lIII Al Therefore, we find that excitation 
of Goldstone modes by a weak rf field will not obscure the 
two-peak signature at zero temperature. At low temper- 
atures of order kT <C ZJ , the initial state still contains 
zero B bosons but does contain a bosons (corresponding 
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to Goldstonc modes) with a bosonic thermal distribution. 
Taking these thermally-excited a bosons into considera- 
tion in the rf signal, we find that I'^^-' is still a delta peak, 
while I^^^ becomes a finite function of \uj — lub\ which 
does not obscure the peak I^^^ . For higher temperatures 
of order ZJ (when the average number of a bosons per 
site is of order 1 or more) , the higher order terms in Hrf 
cannot be ignored. 

APPENDIX C: TIME-EVOLUTION OF A 
SIX-BOSON FOUR-SITE SYSTEM 

The dynamics of finite bosons in a small number of 
lattice sites can be numerically solved without using the 
mean-field approach, and the results reflect the physical 
signatures predicted in our theoretical analysis. Here we 
consider a system of six bosons on four sites located on 
a ring. The ground state is given by 

^(|1)|1) |2)|2) + |1)|2) |2)|1) + |2) |2)|1)|1) 

+ |2)|l)|l)|2)) + i(|l)|2)|l)|2) + |2)|l)|2)|l)) 
+0{J/U), (CI) 



where \ni) |n2) {n^) \n4) means the state with the ith site 
occupied by bosons. The ground state has nonzero 
single-site number fluctuations, which differentiate the 
condensate from the Mott-insulator. Thus the system 
is expected to have an rf spectrum similar to that of 
a many-particle condensate. We use the Bose-Hubbard 
Hamiltonian, plus time-dependent terms which repre- 
sent the interaction with the rf field, to numerically time 
evolve the ground state, and compute the average num- 
ber of bosons transferred from the a state to the b state 
(Fig. [5]). The result shows two peaks with shifts of or- 
der J, as we obtained in Sec. |TT] through the mean-field 
approach. 

To show matter-wave interference, we calculate the 
density contribution per particle of a two-particle four- 
site system (condensate) and that of a four-particle four- 
site system (the Mott-insulator) in free expansion(Fig.[5]). 
Both cases have the same initial density profile but have 
very different interference patterns in expansion. The 
condensate shows more interference peaks and has a 
much higher density in the center at long times than 
the Mott-insulator does. This difference is similar to the 
difference between the phase-coherent state (superfiuid) 
and the phase-incoherent state (the Mott-insulator) in 
the experiment of Ref. Q . 
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